Vortex lattices in the lowest Landau level for confined Bose-Einstein condensates 
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We present the results of numerical calculations of the groundstates of weakly-interacting Bose- 
Einstein condensates containing large numbers of vortices. Our calculations show that these ground- 
states appear to be close to uniform triangular vortex lattices. However, slight deviations from a 
uniform triangular lattice have dramatic consequences on the overall particle distribution. In partic- 
ular, we demonstrate that the overall particle distribution averaged on a lengthscale large compared 
to the vortex lattice constant is well approximated by a Thomas-Fermi profile. 
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I. INTRODUCTION 

Experiments on rotating Bose condensates at high an- 
gular momentum have reached the limit in which the 
vortex cores overlap strongly ^Ij. In this limit, the single 
particle states are restricted to states in the lowest Lan- 
dau level (LLL) 0,13. The vortices cannot be considered 
to interact by pairwise interactions, but have intrinsically 
multi- vortex interactions Q , leading to a very interesting 
regime of vortex physics. 

In an influential theoretical paper studying the nature 
of the vortex lattices in this regime the assumption 
was made that the vortices will form a uniform triangu- 
lar lattice. Under this assumption, it was shown that the 
particle density averaged over a lengthscale large com- 
pared to the intervortex spacing has a Gaussian profile. 
One can, however, expect that in the true LLL ground- 
state the vortices will not form an ideal triangular lattice, 
and that small changes in vortex position could lead to a 
rather different overall density profile; within simple ap- 
proximate considerations, when there are a large number 
of vortices one would expect these deviations to lead to 
a Thomas-Fermi profile for the average particle distribu- 
tion 0, 0, Q , which for a harmonic trap is an inverted 
parabola as a function of the distance from the rotation 
axis. 

We expect that the LLL approximation is an excel- 
lent guide to the low-energy properties (including the 
ground state density profile) whenever the ratio. A, of 
the interaction energy scale to the level spacing of the 
transverse harmonic confinement in an axially symmet- 
ric trap, A = 4:iTh'^asn/(Mhuj±), is much smaller than 
unity (n is the typical number density of bosons, the 
s-wave scattering length, M the particle mass, and uj± 
is the trapping frequency perpendicular to the rotation 
axis). We term A <C 1 the LLL regime. Note that the 
condition for being in this regime does not involve the 
system size. When many vortices are present, it is equiv- 
alent to the healing length being larger than the vortex 
spacing. There will be perturbative corrections to the 



LLL approximation, in powers of the interaction param- 
eter A. An earlier work 0,0] , which extends into the LLL 
regime, finds within a variational ansatz that even if one 
neglects the deformations of the vortex lattice that we 
study here, a TF profile emerges through weak Landau 
level mixing provided the healing length is small com- 
pared to the sample size. Also Ref. jlOj provides a treat- 
ment of the averaged density profile in the vortex lattice 
in the opposite limit in which the healing length is small 
compared to the vortex spacing, and deviations of the 
vortex lattice from uniform triangular were found there. 

In this paper, we provide explicit numerical calcula- 
tions of the groundstates of large numbers of vortices in 
atomic Bose condensates in the LLL limit. By comparing 
our full variational results for the many-vortex ground- 
state with the results obtained under the assumption of a 
uniform triangular vortex lattice, we show that the small 
changes in the vortex position in the true groundstate 
away from a uniform lattice have a dramatic effect on 
the overall density profile. In particular, we demonstrate 
that the coarse-grained particle density of a system con- 
taining a large number of vortices is well approximated by 
a Thomas-Fermi profile. The accuracy of the Thomas- 
Fermi profile in describing this regime can already be 
seen from the results of the numerical studies of Sinova 
et al. where the coarse-grained particle density for a 
small vortex lattice array in the LLL was found to fit an 
inverted parabola. 



II. MODEL 

We consider the rotating Bose gas when interactions 
are sufficiently weak to enter the LLL regime (A ^ 
1) and, furthermore, to allow restriction to the two- 
dimensional (2D) limit in which all particles occupy the 
lowest harmonic oscillator state of the axial confinement 
(accurate provided Airh'^asn/ {Mhu;\\) <C 1, where l^h is 
the trapping frequency parallel to the rotation axis). Al- 
though these conditions are far from being satisfied in 
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atomic gases in the absence of rotation, at high angu- 
lar momentum the radial expansion of the cloud due to 
the centrifugal forces reduces the mean particle density 
and leads the system towards the weakly-interacting limit 
[5, 9]. Indeed," recent experiments [l| have reached con- 
ditions that are within the LLL regime and close to the 
2D limit. 

In this paper we shall restrict attention to the case 
where the filling fraction (the ratio of the number of par- 
ticles to the number of vortices j^J) is sufficiently large 
that the groundstate is a vortex lattice [iTl IT^ and can 
be well-described by Gross-Pitaevskii theory. In the case 
of interest - weak contact interactions - Gross-Pitaevksii 
theory for the groundstate amounts to minimizing the 
interaction energy 
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within the space of states in the lowest subband of the 
axial confinement (z) and in the lowest Landau level of 
the transverse {x—y) motion. Such states may be written 



where ry = {x+iy), and a||_j_ 
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lengths parallel and perpendicular to the rotation axis. 
We choose the normalization constant A such that 
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(3) 



which is why the number of particles N enters in the 
expression for the interaction energy, Eqn.l^. The vari- 
ables rjj are the complex numbers representing the po- 
sitions of the vortices in the x — y plane. Hence, the 
wave function (including therefore the particle density) 
is fully specified by the positions of the vortices. The 
task is to find the positions of the vortices that minimise 
the interaction integral as a function of the total angular 
momentum per particle in units of ?i. 
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Note that, within the LLL states, the combined kinetic 
and potential energy (relative to the zero point energy) 
is proportional to the angular momentum, which is why 
these energies need not be included explicitly in the min- 
imization when angular momentum is constrained. 

The connection between vortex positions and parti- 
cle density lies at the heart of the difficulty of an exact 
derivation of the density profile. Variations in the coarse- 
grained 2D density profile, n2d(a;, y), (the particle density 
|V'(t")P averaged over lengthscales large compared to the 
vortex lattice constant and integrated over the z direc- 
tion) are by necessity tied to variations in the vortex po- 
sitions. (A similar relation of superfluid density [which is 



the same as particle density in the Gross-Pitaevskii the- 
ory] and vortex density as in the LLL approximation 
was found in Ref. 0| .) If one makes the assumption that 
the vortices are only slightly perturbed from a triangular 
lattice, and ignores the variation of the lattice geome- 
try in the energetics, then one is led to the expectation 
that the coarse-grained density profile will be an inverted 
parabola - i. e. Thomas- Fermi like 8] . By imposing the 
normalisation condition (PJ, and noting that for states 
in the lowest Landau level, the angular momentum per 
particle is directly related to the density profile via 
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we find that the Thomas-Fermi profile may be written 



n2d(p) 
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where p = \J x^ + y'^- Thus, the radius of the cloud, R, 
increases with angular momentum per particle as 



R = ^3(£+ 1) ax_ 



(7) 



One can also relate the radius to the angular rotation 
frequency f2 at which the state with angular momentum 
I is stable. Following the approach of Ref. '3| one finds 
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where b = 1.1596 is the Abrikosov parameter for the 
energy of a uniform triangular vortex lattice [l^ and 
.g2D = 2^/27^iVn^as/(a||M). 

On the other hand, if one assumes that the vortices 
form a uniform triangular lattice, then the coarse-grained 
2D density profile is of the form 
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where ny is the areal density of vortices. (For a trian- 
gular lattice of lattice constant a, the vortex density is 
nv = 2/(\/3a^).) Using lO, one finds 
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For the full variational study, a numerical conjugate 
gradient method is used to minimise the interaction en- 
ergy (^) with respect to the vortex positions in ||2Jl and 
subject to the normalisation (j^J for a range of values 
of the fixed angular momentum per particle To this 
end, it is helpful to expand the wavefunction in terms 
of the single-particle states in the LLL with angular mo- 
mentum quantum number m, and write the (normalised) 
wavefunction in the x — y plane as 



m>0 



(12) 



Constraints on angular momentum per particle and nor- 
malisation are then imposed by Lagrange multipliers (in 
terms of physical parameters these are h,{u!± — n)N and 
fiN, where fi is the chemical potential relative to the zero- 
point energy). In this way, the variational equations can 
be expressed in terms of coupled equations for the com- 
plex coefficients Cm- 

As an aside, we note that the resulting variational 
equations can be expressed in real-space form as 

+h{uj^~n) /dV'5(^,^')^0('7') = (13) 

where 

7ra]_ 

is (the integral kernel of) the projection operator to the 
LLL [l|. Notice that \S{r],r]')\^ = e"!''-'''!"/"! /(7rai)2 
falls rapidly for |t? — t^'I > a±. This real-space form of 
the lowest Landau level variational equations is a useful 
starting point for further analytic approximations. 

Here we focus on the numerical solution to the prob- 
lem. In our numerical procedure we work in terms of the 
coefficients c™, which we allow to be non-zero up to some 
maximum angular momentum (< 300) which we check is 
sufficiently large to have no significant effect on the re- 
sults (this is equivalent to setting a maximum number 
of vortices in This procedure is identical to that 

followed in Ref . jj| ; the difference is that here we are in- 
terested in much larger values of the angular momentum 
per particle than the cases studied in Ref. j^l- (For the 
Thomas-Fermi profile, at large angular momentum, one 
expects the number of visible vortices, iVy, to increase 
with the angular momentum per particle as Ny = 3^.) 

III. RESULTS 
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FIG. 1: The particle density for two vortex lattices which 
are calculated as the minimum of the interaction energy Q. 
In the upper entry we have angular momentum £ — 31 and in 
the lower entry £ = 91. The grey-scale code is set to black for 
vanishing particle density and to white for maximum density 
which occurs at the central region of the trap. 



The results of our full variational calculation always 
show a (slightly distorted) vortex lattice. We focus at- 
tention on the case of large angular momentum where 
the lattice spacing appears to be roughly uniform in the 
central region of the condensate. For the range of angular 
momentum that we study, we find that there are stable 
low-energy configurations of the system in which one vor- 
tex is close to the centre of the trap, and around which 
there is an approximate 6-fold rotation symmetry. One 
can then choose one of the lines of reflection symmetry 
to be the x axis. We have tested that such configurations 
are stable to small deviations around the minimising val- 
ues, and that they give the lowest energy that we could 
obtain. 

In Fig. n we show shadow plots of the particle den- 
sity for two examples of the results of the full variational 
study, at ^ = 31 and i = 91. Even though the vortex 



positions of the full variational groundstate appear to be 
very close to the positions of a triangular lattice, there 
are slight deviations from this regular arrangement, espe- 
cially close to the edge of the cloud. These deviations are 
sufficient to give a very different overall density profile to 
that expected from a uniform vortex lattice. 

To illustrate this, we have constructed the wavefunc- 
tions for the ansatz in which the vortices are at the sites 
of a uniform triangular lattice |^. To most closely re- 
produce our full variational results we choose to position 
the lattice such that there is a lattice site at the cen- 
tre of the condensate and there is reflection symmetry in 
the X-axis. The only remaining freedom in the ansatz is 
then the value of the lattice constant, which controls the 
angular momentum per particle. 

In Fig. 121 we show the particle densities for this trian- 
gular lattice ansatz, with lattice constant chosen to give 
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FIG. 2: The particle density for two wavefunctions obtained 
under the assumpion of a perfect triangular vortex lattice. In 
the upper entry we use a lattice constant a = 1.9351ax which 
corresponds to angular momentum £ — 31 and in the lower 
entry we have a = 1.9151ax which gives £ — 91. The grey- 
scale code is set to black for vanishing particle density and to 
white for maximum density which occurs at the central region 
of the trap. A qualitative comparison with Fig. is easy but 
note that the white colour in Fig. does not correspond to 
the same value in the present figure. 

angular momenta £ = 31 and ^ = 91. The differences be- 
tween the particle densities obtained from the full varia- 
tional calculation and from the triangular lattice ansatz 
are large enough to be seen directly from the shadow 
plots. The particle density for the triangular ansatz is 
more peaked at the centre and has a somewhat longer 
tail at large distances. 

To make a more quantitative comparison of the above 
lattices we have determined the angular-averaged parti- 
cle density as a function of the radial distance p from the 
trap centre, by writing rj = pe"^ and averaging the par- 
ticle density over the azimuthal variable tp. The results 
are presented in Figs. |3Ia) and (b). The difference be- 
tween the full variational result and the triangular lattice 
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FIG. 3; The angular- averaged particle density as a function 
of the radial co-ordinate for (a) ^ = 31 and (b) £ = 91. (The 
particle density is in units of a~[^ and p in units of a±.) In 
both figures, the solid line corresponds to the full variational 
result, and the dashed line corresponds to the triangular lat- 
tice ansatz. The dotted lines are the Thomas-Fermi profiles, 
Eqn. 0, for (a) ^ = 31 and (b) £ = 91, while the dot-dashed 
curves are the Gaussian profiles © for the lattice constant 
(a) a = 1.9351ax, and (b) a = 1.9151ax- 



is very clear for i — 31 as well as for ^ = 91. For both 
values of angular momentum, the average density profile 
for the full variational result is very well reproduced by 
the Thomas-Fermi profile, Eqn. © (ignoring the rapid 
oscillations on the lengthscale of the intervortex spacing) . 
On the other hand, the average particle density for the 
triangular lattice ansatz is well reproduced by the Gaus- 
sian formula consistent with the result of Ref. j^. 
The Gaussian and Thomas-Fermi profiles are sufficiently 
different that it can be clearly seen that the Thomas- 
Fermi profile provides a significantly better fit to the full 
variational groundstate than does the Gaussian. 
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FIG. 4: The interaction energy as a function of the 
angular momentum per particle, for the full variational 
groundstate (solid line) and for the triangular lattice ansatz 
(dashed line). [The interaction energy is plotted in units of 

Finally, in Fig. 0] we compare the interaction ener- 
gies per particle for the full variational wavefunction and 
the triangular lattice cases as a function of angular mo- 
mentum per particle. As is required for consistency, 
the full variational wavefunction has an interaction en- 
ergy lower than that of the triangular lattice ansatz (by 
about 10%). The interaction energy for the full vari- 
ational wavefunction is very closely given by the ex- 
pression Ey^^ = 2bg2DN/[9TTa\{i + 1)], which can be 
found by use of the Thomas-Fermi analysis of Ref. I^; 



the corresponding energy for the triangular lattice is 
Ell^l = bg2DN/[4:na\{i + 1)]. This energy reduction by 
a factor of 8/9 at fixed angular momentum, £, is equiva- 
lent to the energy reduction by a factor of 2\/2/3 at fixed 
angular frequency found in Ref. 8]. 



IV. CONCLUSIONS 

In conclusion, we have presented the results of numeri- 
cal calculations of the finite-size vortex lattices in atomic 
Bose condensates at large values of the angular momen- 
tum. We have studied the weakly-interacting limit where 
vortex cores overlap strongly, and the particles are re- 
stricted to states in the lowest Landau level. Our results 
show clearly that the density distribution of particles av- 
eraged over a lengthscale larger than the vortex lattice is 
very accurately given by a Thomas-Fermi profile. These 
results indicate that, in the lowest Landau level limit, 
a fit of the overall density profile found in experiment 
to a Thomas-Fermi profile, Eqn.®, can be used as a 
measure of the angular momentum per particle. Images 
of the groundstate density profile and vortex locations 
can be obtained reliably from the expanded cloud (fol- 
lowing release of the trapping potential), since for states 
in the lowest Landau level the wavefunction of the final 
(expanded) state is directly related to the initial (unex- 
panded) up to a rescaling and rotation . 

While preparing this work for publication, we became 
aware of related numerical studies where similar conclu- 
sions were reached 
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